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FLUX  DIFFUSION  AND  ACCELERATION  IN  A  HIGH  CURRENT 
MODIFIED  BETATRON 

I.  Introduction 

The  evolution  of  a  high  current  electron  beam  in  a  modified  betatron  [1]- 
[4]  during  self  flux  diffusion  and  acceleration  is  described  in  this  paper. 

An  investigation  is  also  made  of  the  effect  on  beam  equilibrium  of  the 
presence  of  a  discrete  set  of  coils  designed  to  compensate  for  flux 
diffusion.  Our  results  are  based  on  calculations  of  equilibrium  solutions  of 
the  cold  fluid  equations  and  their  adiabatic  evolution  as  described  in  earlier 
works  [5,  6].  There  it  is  shown  that  the  fluid  canonical  angular  momentum, 

P0,  characterizes  the  equilibria  because  many  important  fluid  quantities  are 
functions  of  P.,  and  because  surfaces  of  constant  P.  are  drift  surfaces  of 

9  0 

single  electrons.  It  follows  that  beam  equilibria  can  only  exist  on  those  Pa 

9 

surfaces  which  are  closed  within  the  confining  conducting  wall. 

In  Fig.  1  we  show  a  schematic  diagram  of  the  modified  betatron  together 
with  the  experimental  parameters  we  use  (unless  otherwise  specified).  The 
motion  of  electrons  is  broken  into  four  time  scales:  the  poloidal  (r  -  z 
plane)  cyclotron  motion  in  the  toroidal  (BQ)  field  with  period  of  about  10”*^ 
sec,  the  toroidal  major  rotation  in  the  vertical  field  of  about  10~®  sec,  the 
poloidal  drift,  vp,  in  the  B0  field  of  about  10“  7  sec,  and  finally  the  time 
scales  for  acceleration  and  flux  diffusion  on  the  order  of  10“^  -  10“^  sec. 

In  contrast  to  concurrent  studies  [7]-[10]  using  particle  simulations,  our 
fluid  description  averages  out  the  individual  motion  of  electrons  and  studies 
the  adiabatic  development  of  the  beam  on  time  scales  longer  than  the  poloidal 
drift  period. 

In  our  axlsymmetrlc  system  at  steady  state,  all  the  forces  on  the  beam 
are  in  the  poloidal  plane.  The  outward  Y?0/r  centrifugal  force  is  balanced  by 
the  Inward  v0  x  Bz  force  and  the  defocuslng  effects  of  the  self  electric  field 
and  image  charge  are  balanced  by  the  focusing  forces  due  to  the  external 
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Index  ri,  the  poloidal  self  magnetic  field  with  image  currents  (v  x  B  ),  and 

o  P 

the  v  x  B.  force.  In  a  high  current  betatron,  the  repulsive  electrostatic 
P  0 

forces  dominate  the  attractive  forces  due  to  external  focusing  and  magnetic 

focusing.  Thus,  the  additional  v  x  B  is  necessary  for  equilibrium  and  the 

P  “ 

resulting  electron  drift  is  diamagnetic.  The  other  forces  in  the  poloidal 
plane,  the  inertia  and  centrifugal  forces  due  to  poloidal  motion  are  neglected 
in  our  formalism  because  they  are  negligible  when  the  toroidal  field  is 
sufficiently  large  and  the  density  and  energy  sufficiently  small.  These 
effects  will  be  discussed  in  a  future  work. 

Since  the  system  is  assumed  to  be  axisymmetric,  P  is  conserved,  and  the 

0 

key  to  calculating  the  equilibrium  is  the  calculation  of  the  P.  surfaces.  It 

2 

can  be  shown  that  the  fluid  energy  E  «  ymc  +  q<j>  and  g  =  r  BQ  are  constant  on 
the  PQ  level  curves,  and  that  the  Pg(r,z)  surfaces  can  be  found  by  specifying 
the  dependence  of  E  and  g  on  Pa .  For  parameters  which  vary  slowly  compared  to 
the  poloidal  drift  period,  the  evolution  of  the  beam  may  be  considered 

adiabatic.  The  adiabatic  constants  of  the  motion  are  P.  itself,  and  the 

0 

toroidal  flux  enclosed  by  a  Pa  curve,  $,.(PQ).  These  two  constants  of  the 

0  n  0 

motion  determine  the  evolution  of  E(P  ,t)  and  g(PQ ,t) ,  and  thus  generate  an 
evolving  series  of  Pg  surfaces. 

In  a  previous  work  [5,6],  the  behavior  of  the  equilibria  during  flux 
diffusion  and  acceleration  was  presented  for  a  modified  betatron  with 
rectangular  boundaries.  In  this  paper  we  present  results  for  a  circular 
boundary  using  a  modified  version  of  the  original  iterative  computer  code  in 
which  the  outer  loop  iterations  have  been  accelerated  with  a  fast  direct 
elliptic  solver.  We  find  no  radically  different  results  in  this  new 
geometry.  Two  general  conclusions  we  reach  are  that  the  effect  of  flux 

diffusion  is  much  less  disturbing  to  the  ?a  topology  in  the  circular  geometry, 

0 


and  that  the  window  of  acceptable  field  indices  for  acceleration  is  wider 
because  of  the  shape  of  the  new  confining  wall. 

After  injection  and  trapping,  the  beam  in  the  modified  betatron  is 
accelerated  in  three  phases:  (1)  preacceleration,  (2)  diffusion,  and  (3)  main 
acceleration.  We  will  deal  with  the  diffusion  and  acceleration  phases  since 
they  occur  on  time  scales  of  at  least  two  orders  of  magnitude  greater  than  the 
poloidal  drift  period.  In  our  discussion  of  diffusion,  we  will  examine  the 
effect  of  Imperfect  flux  diffusion  compensation  on  beam  equilibria.  A 
discussion  of  the  numerical  solution  procedure  and  the  new  elliptic  solver  can 
be  found  in  the  appendix. 

II.  Flux  Diffusion 

a)  Constant  Energy 

Flux  diffusion  is  simulated  in  our  system  by  specifying  the  boundary 
conditions  on  the  flux,  -  rAQ,  as  follows: 

♦(r,s)  -  <pe(r,z)  +  e  *r(r,z)  +  i|>o  (1) 

where  is  the  well  known  flux  due  to  a  thin  ring  carrying  beam  current  I 

[11],  is  due  to  the  external  vertical  field  coils  and  the  constant  <po  is 

determined  so  that  P.  at  the  reference  orbit  is  conserved.  The  constant 
e  reflects  the  degree  of  diffusion;  e  ■  0  represents  no  flux  diffused  and 
e  ■  l  full  flux  diffusion.  The  flux  \|)Q  is  the  additional  flux  through  the 
major  orbit  that  must  be  provided  by  external  coils.  In  this  set  of 
calculations,  the  energy  of  the  beam  is  held  constant  while  the  vertical  field 
is  adjusted  appropriately  to  hold  the  beam  at  the  initial  equilibrium  radius 

rQ.  In  the  first  of  these  runs,  rQ  ■  100  cm,  the  field  index,  t\  -  0.5,  and  as 
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shown  in  Fig.  2,  we  find  that  the  vertical  field  Bv  increases  29%  from 
Bv  ■  120  gauss  with  no  flux  diffused,  to  By  *  155  gauss  with  all  the  flux 
diffused;  »|>  also  increases  during  diffusion  from  0  to  2.36x10*'  gauss-cnr,  but 
ip(ro)  remains  almost  constant.  We  shall  explain  the  Interdependence  of  these 
variables  shortly.  In  Fig.  3  we  show  the  P.  surfaces  before  and  after  flux 
diffusion.  The  x-point  at  approximately  90  cm  hardly  moves  toward  the  beam  at 
all,  so  there  is  no  threat  of  transition  to  the  paramagnetic  curves  outside  of 
the  separatrlx.  This  then  contrasts  to  our  findings  with  rectangular  geometry 
where  the  separatrlx  moves  in  during  diffusion,  ultimately  disrupting  the 
equilibrium. 

If  flux  diffusion  were  to  occur  with  a  centered  beam  and  constant 
external  fields,  the  beam  would  lose  energy  and  its  equilibrium  radius  would 
decrease  because  the  decrease  in  equilibrium  radius  associated  with  the 
reduction  of  yq  Is  greater  than  the  increase  in  equilibrium  radius  associated 
with  the  enhanced  hoop  forces  [7].  In  our  calculations,  energy  and 
equilibrium  radius  are  fixed,  so  the  vertical  field  Increases  to  prevent  the 
increase  in  rQ  produced  by  the  enhanced  outward  hoop  and  image  forces. 

In  the  next  run  we  set  rQ  -  95  cm  and  find  that  By  decreases  from 
Bv  -  121  gauss  to  Bv  -  81  gauss  during  flux  diffusion,  while  \|>q  increases  from 
0  to  6.62x10^  gauss  -  cm^  (see  solid  lines  in  Fig.  2).  The  Pa  surface  after 
diffusion  is  shown  in  Fig.  4.  This  time  the  separatrlx  does  move  closer  to 
the  beam  from  about  4  cm  away  to  1  cm  away  when  flux  has  diffused,  but  again 
the  beam  is  not  disrupted. 

The  explanation  for  the  decrease  in  Bv  is  connected  to  the  shift  in  the 
self  magnetic  field  during  diffusion.  As  the  repelling  image  currents  in  the 
wall  decay,  the  net  forces  on  the  beam  point  inward  due  to  the  electrostatic 
attraction  to  the  inner  (nearer)  wall.  Thus  the  Imposed  inward  force,  that 


from  Bv,  must  also  decrease  to  maintain  rQ.  The  external  flux  ($o)  again 
increases  to  compensate  for  the  loss  of  self  flux  and  By  flux  at  the  beam. 

Thus  the  question  of  how  By  must  change  during  flux  diffusion  is  extremely 
sensitive  to  beam  position. 

b)  Constant  Vertical  Field 

In  another  set  of  runs  we  diffuse  the  self  flux,  but  fix  By  and  allow 
the  energy  to  adjust  during  diffusion.  As  in  the  fixed  energy  results,  the 
separatrlx  does  not  move  closer  to  the  beam  when  rQ  -  100  cm,  but  comes  within 
0.25  cm  of  the  beam  when  rQ  ■  95  cm.  Since  the  separatrlx  in  the  later  case 
hits  the  wall  of  the  torus,  this  is  a  potentially  dangerous  situation,  and 
beam  disruption  is  likely  to  occur  for  beam  with  rQ  -  95  cm.  In  Fig.  5  we 
show  the  effect  on  beam  energy  yQ,  on  flint  'l»(rQ)  and  on  the  external  flux 
of  flux  diffusion  with  constant  vertical  field  By.  At  rQ  -  100  cm,  the  energy 
decreases  during  diffusion,  from  yq  ■  3.8  to  y  »  2.0,  to  counter  the  Increased 

A 

outward  (in  the  r  direction)  forces  on  the  beam.  At  rQ  -  95  cm,  the  energy 
increases  during  diffusion,  from  y  -  4.3  to  y  -  6.45.  This  time  the  net  force 
from  the  image  currents  is  outward  and  its  loss  must  be  compensated  by 
increased  centrifugal  forces. 

When  the  flux  diffuses  at  constant  energy,  \f>(rQ)  remains  nearly  constant 
but  By  changes  dramatically  (Fig.  2).  When  flux  diffusion  at  constant 
By,  yQ,  and  t(r0)  both  change  dramatically  (Fig.  5).  In  fact,  during  flux 
diffusion,  our  results  show  that  yQ  and  4*(ro)  are  related  linearly,  exactly  as 
they  are  during  acceleration.  This  follows  because  Pa  •  rmyv  -  e^i  is  a 
constant  of  the  motion,  and  r  and  v.("  c)  are  almost  constant  during  flux 
diffusion.  The  behavior  of  i|io  (which  must  be  supplied  externally)  in  Figs.  2 
and  5  can  be  directly  Inferred  from  the  behavior  of  $(rQ).  The  change  in  flux 
at  beam  center  is  produced  by  four  effects 
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A*(r0)  *  Wv  +  Wf  +  Wc  +  A<|/q  (2) 

whera  ^  is  the  flux  due  to  the  vertical  field,  ^  is  due  to  self-flux 
diffusion,  is  due  to  the  beam's  current  density,  and  t|>o  is  the  constant 
external  flux.  Since  our  current  densities  remain  nearly  constant,  Ai|»c  will 
be  neglected  in  our  simple  argument.  During  flux  diffusion  at  constant 
energy,  i|i(ro)  is  almost  constant  and  ij>o  maintains  \p(ro>  by  compensating  for 
the  loss  of  self-flux  and  the  loss  or  gain  of  flux  from  the  increasing  or 
decreasing  Bv  field  (Fig.  2).  During  flux  diffusion  with  constant 
Bv,  i|»(rQ)  is  tied  to  yQ  and  again  must  supply  the  necessary  flux  according 
to  Eq.  (2)  (Fig.  5). 

III.  Flux  Diffusion  Compensation 

Because  it  may  not  be  possible  to  adjust  the  vertical  field  fast  enough 
to  maintain  equilibrium,  and  furthermore,  because  the  proper  Increase  in  By 
or  ipQ  with  e  is  extremely  sensitive  to  the  beam  position  in  the  liner, 
it  has  been  proposed  that  a  discrete  set  of  coils  be  used  to  compensate  for 
the  decay  of  eddy  currents.  In  this  section  we  simulate  the  results  of  an 
error  in  this  process  of  compensation. 

He  suppose  that  the  set  of  external  conductors  is  placed  along  the  minor 
cross  section  of  the  torus  having  a  pololdal  distribution  resembling  the 
distribution  of  wall  currents  in  a  perfect  conductor.  We  suppose  further  that 
the  error  made  in  this  compensation  produces  a  sinusoidal  variation  of  flux  at 
the  boundary  (the  flux  at  the  boundary  would  be  identically  zero  if 
compensation  were  perfect).  Then  the  boundary  condition  on  t  becomes 
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#  «  *e  +  s  *  +  #q  (3) 

where  t  oscillates  along  the  walls, 

#(r,z)  -  #  cosl6# ,  (4) 

m 

i|»m  is  the  maximum  value  of  #r  and  #  signifies  poloidal  angle. 

Figure  6  shows  the  results  for  e  -  0.25  and  1.0.  These  curves  should  be 
compared  to  Fig.  3a  which  corresponds  to  e  ■  0.  The  presence  of  the  coils  is 
markedly  seen  in  Fig.  6a,  but  the  distance  between  the  separatrix  and  the  beam 

i 

decreases  only  by  1.5  cm  between  e  ■  0.25  and  e  ■  1.0,  and  the  separatrix 
remains  well  away  from  the  beam.  Furthermore,  and  in  contrast  to  our  previous 
constant  energy  flux  diffusion  results,  the  vertical  field  Bv  increases  only 
0.83Z  between  e  ■  0  and  e  ■  1.  These  results  are  very  optimistic  in  that  they 

! 

show  that  the  beam  is  not  affected  by  flux  diffusion  as  long  as  the  diffusion 
is  compensated  by  the  set  of  discrete  coils  mentioned  above.  Furthermore,  the 
beam  is  not  sensitive  to  the  level  of  discretization  error. 

IV.  Poloidal  Error 

Next  we  consider  the  possibility  that  an  error  is  made  in  the  poloidal 
distribution  of  currents  in  the  set  of  coils  described  above.  Suppose  this 

j 

j 

error  is  represented  with  the  following  boundary  condition  on  #  j 

| 

I 

j 

♦(r,z)  ■  1>e  +  ?  +  e  cos#  +  (5) 
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where  again  is  the  maximum  value  of  flux  produced  by  a  thin  ring  of  current 
carrying  beam  current  I,  and  $  signifies  poloidal  angle.  As  shown  in  Fig.  7, 
we  find  that  both  the  vertical  field  necessary  to  maintain  the  equilibrium 
radius  and  the  external  flux  are  very  sensitive  to  these  poloidal  errors, 
and  that  this  sensitivity  increases  with  beam  current.  Equivalently  one  would 
find  that  for  a  fixed  external  field,  the  equilibrium  radius  is  very  sensitive 
to  the  degree  of  poloidal  error.  The  percentage  change  in  Bv  for  5%  poloidal 

error  is  18Z  for  I  ■  10  kA,  8Z  for  I  -  3  kA,  and  3Z  for  1*1  kA.  Thus,  when 

compensating  for  flux  diffusion  with  external  coils,  it  is  essential  that 
poloidal  errors  in  current  distribution  in  these  coils  be  very  small. 

V.  Acceleration 

To  treat  acceleration  of  the  beam  we  assume  perfect  flux  diffusion 
compensation,  and  increment  EQ,  the  energy  at  the  equilibrium  radius  (or  o- 
point  of  the  PQ  surface).  As  in  the  rectangular  geometry  we  find  that  because 
of  the  increased  focusing  forces,  the  separatrix  moves  in  toward  the  beam  as 

its  energy  Increases.  Ultimately  the  separatrix  hits  the  outer  edge  of  the 

beam  and  the  transition  from  diamagnetic  to  paramagnetic  motion  occurs.  The 
beam  will  almost  certainly  not  survive  this  transition  if  the  separatrix 
intersects  the  liner  wall. 

We  now  find  that  as  yq  increases  from  about  4  to  10  (the  transition 
energy),  the  orbit  is  confined  during  the  transition  for  0.4  <  n  <  0.6,  but  is 
not  confined  for  field  index  n  -  0.3  or  ri  ■  0.7.  The  P  topology  for  two 

field  indices,  n  ■  0.5  and  n  ■  0.7  is  shown  in  Fig.  8  just  before  the  edge  of 

the  beams  undergo  transition.  In  the  n  -  0.5  case,  the  energy  yq  *  9.26,  and 
the  separatrix  is  well  contained  in  the  liner.  In  the  q  -  0.7  case,  the 

energy  yq  ■  8,  and  the  separatrix  (in  this  case  a  hidden  level  curve  with  two 
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x  points  at  r  -  99  cm,  z  -  ±  4  cm)  runs  into  the  liner. 

At  this  point  a  few  other  considerations  are  worth  mentioning.  We  have 
run  some  acceleration  calculations  at  beam  current  of  1  kA  and  found  that  the 
beam  will  survive  the  paramagnetic  transition  for  a  broader  range  of  field 
indexes  (almost  as  high  as  0.3  <  n  <  0.7).  We  have  also  investigated  the 
consequences  of  initializing  the  beam  with  a  much  flatter  density  profile  than 
usual  and  found  no  effect. 

VI.  Summary 

The  effects  of  flux  diffusion,  inadequate  flux  diffusion  compensation, 
and  acceleration  on  the  beam  in  the  modified  betatron  have  been  examined  in 
this  paper.  Our  formalism  assumes  that  the  poloidal  drift  in  the  toroidal 
field  is  well  below  the  Brillouin  limit  ("slow  mode"),  that  the  beam  is 
"cold",  and  that  the  parameters  affecting  the  beam  change  adiabatically.  We 
find  that  beam  equilibrium  can  be  lost  during  flux  diffusion  unless  the 
diffusion  is  accurately  compensated.  We  also  find  that  a  10  kA  beam  will 
survive  the  diamagnetic  to  paramagnetic  transition  provided  the  field  index  is 
within  the  limits  0.4  <  n  <  0.6. 
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Fig.  1  Modified  betatron  schematic  with  typical  parameters 


Z-|  =  12.5  cm 
Rj  =  87.5  cm 
R2  =  112.5  cm 
r0  =  100.0  cm 


%  FLUX  DIFFUSED 

Fig.  2  Flux  diffusion  at  constant  energy. 


Mi  x  10s GAUSS- 


-  GAUSS 


Appendix;  Description  of  the  Iterative  Equilibrium  Code 


In  this  appendix  ve  discuss  the  equations  governing  the  self-consistent 
equilibrium  of  a  modified  betatron,  and  give  an  outline  of  the  solution 
procedure.  Our  method  is  reminiscent  of  that  used  by  Grad  for  calculating  the 
adiabatic  compression  of  a  two-dimensional  plasma  [12-14]. 

The  cold  fluid  equation  of  motion  for  the  beam  in  equilibrium  is 


v  •  Vp 


q(v  *  |  —  v*) 


(6) 


where  p  *  ymv,  y 
induction,  and  $ 


2,  2 


1/2 


is  the  relativistic  factor  (1-v/c)  ,B  the  magnetic 
the  scalar  potential.  By  toroidal  symmetry  we  can  write 


B  -  7i>  x  78  +  g70,  .  (7) 

where  i|>  *  rA.  is  the  poloidal  flux,  A.  the  toroidal  component  of  the  vector 
potential,  g  -  r  BQ,  and  0  the  toroidal  angle.  Using  ampere's  law  and  the 
toroidal  component  of  Eq.  6,  we  can  show  that  the  level  curves  of  the  toroidal 
canonical  angular  momentum 


PQ  -  ymrv0  +  q*/c 


(8) 


constitute  drift  surfaces  of  the  electrons  in  the  poloidal  plane  and  that 


g  -  g(p0)* 


(9) 


If  poloidal  inertia  and  centrifugal  force  are  ignored,  the  poloidal 


component  of  Eq.  (6)  Implies  that  the  energy 


Is  a  function  of  PQ,  and  that 


..  8_ 

4imr2 


dg_ 


dP0 


dE 

dP9 


Q<V’ 


(U) 


where  n(r,z)  is  the  density. 

Defining  K(PQ)  ■  g  ^  ,  we  can  rewrite  Eq.  (11)  for  the  density. 


n 


K<V  # 

4irr2 [v0/r  -  Q(P0)] 


'  (12) 


Finally  we  compute  the  electrostatic  potential  $  and  pololdal  flux 
f  according  to 


V24  -  4irne  (13) 

* 

7  -  4irnery0/c,  (14) 

*  2  -2 

where  7  t|>  ■  r  7  •  (r  7^),  and  the  boundary  conditions  for  acceleration  are 
those  of  a  conducting  wall;  <p  •  0  and  »  i|»  +  \|>  where  ^  is  the  flux  due  to 
the  external  vertical  field  coils  and  t|io  is  a  constant.  For  flux  diffusion  we 
model  the  flux  on  the  boundary  to  be 
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where  pf  is  the  well-known  flux  due  to  a  thin  ring  carrying  beam  current  I, 
and  e  Is  a  measure  of  how  much  flux  has  diffused.  For  both  acceleration  and 
flux  diffusion,  pQ  is  an  externally  imposed  constant  flux  chosen  so  that  PQ  at 
the  reference  orbit  is  conserved.  pQ  must  be  supplied  by  an  external  coil  to 
maintain  equilibrium  during  flux  diffusion  or  acceleration,  and  is  defined  to 
be  zero  at  the  first  equilibrium. 

Elliptic  Eqs.  (13)  and  (14)  are  solved  using  a  special  fast  Fourier 
technique  which  is  able  to  handle  circular  (or  more  general)  boundaries.  The 
general  idea  is  to  solve  the  elliptic  equations  in  a  rectangular  region 
enclosing  the  actual  boundary.  In  this  rectangular  region  the  usual  technique 
of  Fourier  transforming  in  one  direction  and  solving  the  resulting  tridiagonal 
system* in  the  other  direction  works.  Next  a  "Green's  function"  matrix  is 
generated  which  computes  the  discrete  surface  charges  on  the  real  boundary 
necessary  to  produce  the  desired  potential  (or  flux)  there  [15-17].  This 
approach  can  also  be  interpreted  as  a  method  of  permuting  the  matrix  obtained 
by  standard  finite  difference  in  the  disc  to  a  matrix  on  which  the  standard 
FIT  methods  will  work  [18].  In  this  elegant  framework,  Neuman  boundary 
conditions  on  Irregular  regions  can  also  be  dealt  with  easily. 

In  our  computer  code,  Eqs.  (8)  -  (14)  are  solved  iteratively  for  the 
given  experimental  parameters  such  as  beam  current,  radius,  energy,  etc.,  and 
for  given  values  of  boundary  conditions  on  4  and  p.  Beam  equilibrium  at  any 
point  later  in  the  acceleration  or  flux  diffusion  cycle  can  be  found  merely  by 
specifying  the  energy  of  the  beam  or  the  boundary  conditions  on  p.  The  new 
equilibrium  is  found  by  conserving  the  adiabatic  constants  of  the  motion. 
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